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ABSTRACT 


The praetical and rigorous solutvons of the poecmenan 
problem associated with the harmonic oscillation of a sub- 
merged rigid body of arbitrary shape is presented. The use 
of Green's function reduces the determination of the poten- 
feral to the solution of an integral equation. The integral 
eguatlon is solved numerically and the dependency of the 
hydrodynamic quantities such as added mass and damping 
MErerlTelcncs Of the object on the frequency of the oseriia-— 
tion is established. 

Pcvctalbeem@cexs are made on the numerical results sees 
include the Haskind's relations check, an energy check for 
the radiation problem, and comparisons with closed form 
solutions. All the checks and comparisons are successful and 
it appears that the numerical procedure employed yields valid 


ema accurate results. 
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dymamieG, presstimenassOCia tedmmaatn the) 


NOMENCLATURE 


Description 
characteristic length of object 
wavelength parameter, a = 2na/L, dimensionless 


A A . Sth 
wave force or moment coefficient in the i 


direction, dimensionless 

elementary dimensional surface area 
elementary dimensionless surface area 
energy transmitted over one period 
diSteibutilon pune Gren diem cmen hess 


° -th 
wave force or moment component in the i 


direction, dimensionless 

acceleration of gravity 

res My sell Ihe tenon 

water depth 

relative depth of water, dimensionless 
functions defined by (19), dimensionless 
wave height 

imaginary part 

unit normal vectors, dimensionless 


dynamic pressure associated with the wave 
interaction with the fixed object 


th 


mode of oscillation, dimensionless 
we JLUisal aie Ikeve ie Wetec 

Teal pame 

plan polar coordinate 


plane polar coordinate, dimensionless 








Desc mil lous 


horizontal distance between points (x,y,Z) 
and (§,n,%t), dimensionless 


surface area of the object 


surface area of the object, dimensionless 


time 


wave period 


velocity potential for incident wave, dimensionless 


velocity potentials for the radiation and scatter 
problems, dimensionless 


amplitude of motion in gee 


Girect lon 


amplitude of linear or angular motion, dimensionless 


Spacial variables 


Spacial variables, dimensionless 


gta component of force or moment coefficient 


associated with the jt 


component of oscillation 


of the object 


qe component of the dynamic force or moment 


due to the jth 


mode of oscillation 


added mass coefficient 


added mass coefficiemt, dimensionless 


damping coefficient 


damping coefficient, dimensionless 


Darae 


delta function 


rectangular Cartesian coordinates, dimensionless 


plane polar coordinate 


angular displacement of oscillating object 


o*a/¢ 


(eee 


2 


= a tanh (ah) 


1/2 
) 








Subscripts 


1,) 


Description 
angular frequency, o = 21/T 
velocity potential for incident wave 


velocity potential for oscilieg wonmof jobject 
OT eSecacter 


total velocity potential for diffraction problem 


used to denote either the mode or oscillation 
(or Scatter) or the dincectiomeo: 4 Troree iF 
moment component, depending on context 
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I. “CENTROBUCT LON 


The need for the utilization of tie natumatenesources 
mcdded in coastal waters has led to ehemrecdster create, 
knowledge of the interaction of surface gravity waves with 
large submerged objects such as oil storage tanks and pipe- 
menes. A new problem arises in the case of large bodies 
immersed in shallow water, the scattering of the incident 
mame and the free surface effect become an important considera- 
mem in the calculation of wave forces on fixed objects and 
hydrodynamic forces produced by oscillation of the body in 
otherwise still water. 

The analysis which considers these two factors has come 
to be known as diffraction neces In this approach separa- 
tion and viscous effects may be neglected, as was confirmed 
by Sarpkaya and Garrison [1], provided the amplitude of the 
foeeton 1S sufficiently small. The problem is set up in 
terms of a velocity potential and the potential function 1s 
Metermined by use of a source distribution. The pressure 
on the immersed surface can be found once the source strengths 
are known by use of Bernoulli's equation and the resulting 
monees and moments are determined from the pressure distribu- 
mom by Surface integration. 

Recently, Garrison and Seetharama Rao [2] worked out an 
analysis for wave interaction with submerged objects based on 


Minear diffraction theory and applied the analysis to calculate 


isl 








wave forces acting on a submerged hemisphere. They showed 
that for small values of the wave height to sphere diameter 
Gatio, the wave force coefficients can deme s temmecscnmed 
by the diffraction theory and correlated as functions of the 
parameters, 2na/L and h/a where h denotes the water depth, a 
the sphere radius, and L the wave length. 

Aside from the wave force problem, an additional problem 
may be considered which is referred to as the "radiation 
problem.'"' This involves the effect of a rigid body in calm 
water undergoing forced harmonic oscillations in its various 
degrees of freedom. Of interest here are the hydrodynamic 
jeeces which may be characterized by the coefficients of 
added mass, added moment of inertia, and damping. Although 
often solved independently, the diffraction and radiation 
problems may be solved concurrently because they are mathe- 
matically similar. 

Relations which eliminate the need to solve the problem 
of wave diffraction past an object by solving the forced- 
oscillation problem in still water were developed by 
ieee. Haskind in 1957. They relate the force and moment 
Semponents of the diffraction problem to the asymptotic 
welocity potentials for the corresponding radiation problem. 
That is, if the waves produced by the oscillation of an 
object in still water are known, Haskind's relations may be 
used to obtain the wave forces (or moments) produced by 


Surface waves interacting with the same object held fixed. 


IBZ 





Conservation of energy dictates that aenalancee snus tecrrse 
between the energy required to oscillate an object, and the 
wave energy transmitted by surface waves across some control 
volume surrounding the object but at a large radial distance. 
The energy required to oscillate the object is directly 
related to the damping coefficient, and the damping coefficient 
May be obtained by integration of the pressure over the immersed 
surface. In view of conservation of energy this same energy 
mepedrs at a large distance from the object in the form of 
manated wave energy. Therefore, the dampine coefficient may 
be obtained by accounting for this energy flux. The compari- 
amor Chis coetfiici1ent obtained trom these two methods, 1.c.. 
the near field and far field solutions, provides a valuable 
check on the accuracy of the numerical results generated. 

The application of Haskind's relations and energy considera- 
mmems tO calculate wave forces on large submerged objects is 
the purpose of this thesis. A computer program was developed 
to carry out the numerical calculations and to compare the 


mesults with those from diffraction and radiation theory. 
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II. THEORETICAL CONSIDERATIONS 


A. A REVIEW OF WAVE FORCES ON SUBMERGEDSOLGEGIs 

Consider a rigid body of arbitrary shape having a charac- 
teristic size a submerged in water of depth h as illustrated 
in Fig. 1. It is assumed that the amplitude of the motion 
of the fluid is sufficiently small in relation to the size, 
a, that viscous effects can be neglected. The object may 
@eemay not be in contact with the free surface or the bottom- 
tne problem set forth herein deals with the fluid motion and 
forces induced by the small amplitude oscillation of the 
object in its six degrees of freedom as well as the fluid 
motion associated with the interaction of the fixed object 
with a train of regular waves. 

The small amplitude oscillatory motion of the body with 
frequency about its equilibrium position is described by the 


metationships 


= ~ Oo -4 
X(t) = X, Re[e TOL a = Is (1a) 


.) 


6. Rel[e 


-iot 
; ] 


L/P) 
me 
ct 
Ww 
I 


5 ES Gy (1b) 


where X. denotes the amplitude of motion in translational 
escillation and 0. denotes the amplitude of motion in rota- 
mional oscillation. The subscript 1 = 1, 2, 3 corresponds To 
oscillation in the x, y, z direction, respectively, and 


1 = 4, 5, 6 corresponds to angular oscillation about axes 
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YdIaxS uot Turzaq 
(2 “yz+u-°3) 


[[T ein3sty 
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passing through the origin of the obj) eC tiem cl etoecic 

X, Y, zZ axeS, respectively. The second problem considered 

here is the interaction of a train of regular surface waves 
with the object fixed in space. the incidént®waves of wave 
height H and wave length L are assumed to progress in the 
positive x-direction and the fluid motion is assumed to be 
incompressible, irrotational, and harmonic with time dependence 
oo imeall “Gases, @Gonsequently sa Vvelocityepotengialme tse 


such that the fluid velocity vector may be defined as 


qs - Re[Vo,e 77°} ee as (2 


where d= denotes the velocity potential associated with the 
Homerton induced by oscillations in the six degrees of freedom, 
In the case of regular wave interaction with the fixed object, 


the velocity potential is given by the sum 
b1 = 6 + 6, (3) 


where 5 denotes the velocity potential of the incident wave 
in the absence of the apHIEcE and 4 denotes the velocity 
potential of the scattered wave due to the presence of the 
rigid body. For this case, the fluid velocity vector is 


given by 
q' = Rel[¥(o, + oe OT]. (4) 


Having assumed irrotationality of the fluid, >; must satisfy 
the continuity equation which takes the form of the Laplace 


equation, Wore 0 (5) 
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within the. flue cequen, 
Assuming that the amplitude of the motion is small, the 
velocity squared terms in Bernoulli's equations may be neg- 


lected. Therefore, the dynamic pressure is given as 


P, = Re[oop,e 7"), 5 = 1, 2,...6 (6) 


For the second problem, involving wave interaction with the 


fixed object, the pressure is given by 


Ppt = Relea (>. it be 77] : (7) 


The velocity potential must satisfy certain boundary 
conditions in addition to Eq. (5). These include the 


linearized free surface boundary condition, 


od. 2 
Bern Oy cz) = — ¢. (x, 0, Ze jay 2 eae (8) 
ee g 4 


as well as the kinematic boundary condition on the bottom, 


3. a 
—2, (&.-h, 2) = 0 , i= 1, 2,23 (9) 


ay 
Since we satisfies both Eqs. (8) and (9), it follows, there- 
fore, that b4 must also satisfy these conditions. 
In addition, 5 must satisfy the boundary condition on 
the surface of the body defined by s(x, y, z) = 0. This con- 
dition states that for an impermeable boundary, the normal 


velocity of the fluid must just equal the velocity of the 


Burkace normal to itself. In the case of oscillation OCfVvtne 








rigid body, the kinematic boundary conditions on the surface 


of the object take the form 


9>4 °o 
- (x, Was Z y te) = Xin, 
on 
ace yk 
—— (x, y, z, t) = X5n 
on 2°y 
9b ve 
ares Ce ee) al, 
on 
(10) 
2%, _—— _ 
eas ee 1) O,tth + y)n, S xny | 
nN 
9¢¢ re) = ms 
= (on OG Be 18) 8-[zn, = xn, |] 
on 
996 e a _ 7 
se Cos 85 AG) Tol ay Sn y)n,] 


ere n = in, + jny + kn, denotes the unit normal vector on 

tme surface of the object directed outward into the fluid. For 
the case of the interaction of regular waves with a fixed, 
rigid body the normal velocity must be zero which may be 


written 


1G eee 20 eT ae 
=== UE: ne a —— Cr eee ee (ely) 
on on 

Finally, the velocity potential must satisfy the radiation 


condition which permits only an outgoing progressive wave at 


a large distance. Hence it is required that o. have the 
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asymptotic form which satisfies 


Sa) 2 Pia -_ 
.(F, 0, 7) - AC),  SOshlkUhey)) qt(kry-ot) , 9 as z+ @ (12) 
J J cosh(kh) 


joe- 1, 2, «2:7 where tT) - (x? + z2yll2 and = sare tana) oo) 
where B,(8) is some unknown complex function of 9G. 

The wave number is defined as k = 21n/L where L is the wave 
length and is related to the frequency of disturbance according 


to 
os : 
ro kK tanh (kh) . (13) 


The velocity potential of the incident wave alone progres- 


feao ein the positive x-direction 1S given by 


— Te igi cosh[k(h + y)] ,-ikk e 
cosh (kh) 

where 5 = H/2 represents the amplitude of the incident wave, 

H being the wave height. 

In carrying out the solution and to show more clearly the 
dependence of the solution on the relative wave length parame- 
ter, a = 2na/L = ka, and relative depth, h = h/a, a being a 
typical length of the body, it is convenient to make the space 


variables and amplitudes dimensionless, i.e., 


SS eR oe SR 8 we a= ea Se Saye 


1) 





and then to introduce the  onecsure stimetaon be by 


iog; (x, y, z)/gax, =| i iegncda (aun) u, (Xig-0¥ 5 ze = ee eee 


ino (x, Was z)/gan -aul (x, Y> zy) 








(15) 
iodo (x, y, z)/gan = -auy(x, y, Zz) 
The complex dimensionless dynamic pressure amplitude can now 
be written by use of the linearized form of Bernoulli's 
equation as 
P; = a tanh (ah) u, (Xs Vo. Zest le eee 
A Ce) 
oe = cosh[a(h + y)] e er -au, (x, Vie Z) 
cosh(ah) 
where the complex amplitude of the pressure p is defined as 
Ip ; 
1, = Re[p, e *°°] 
ogaXx. J 
J 
(17 
7 
= = Refpe 2") 
Ogan 
Tae boundary value problem which descripes tne fluid 
Hoeraon arising from the oscillation of the rigid body in its 
Six degrees of freedom as well as the scattering of the 
mierdent wave can now be written concisely in terms of dimen- 
Suenless parameters. The potential us (x, V6 92) 5 7 = iL3 2a 
continuous in the fluid region is sought such that 
Veu, (x, ¥, zZ2) = On ye <8 (ayy Ges 
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syr (x, 0, 2) - a tanh (ah) u, (x, 0, z) = 0 (b) 
alls 
Sa (x5 =h, 2) = 0. outside ls (ee (Ee) 
ou. 
gal OX ys 2) = hy (x, y, 2) on s(x, y, 2) = 0 (A) 


=1/2 Coshjath. = y )| seetane 


Meer, 8, y) - bee) vr Sahel = i mas Tote 


where Yale) Ps Some dimens Vonless “comp lex cire & loimo ie 
meres. zz) = 0 represents the surface of the object 1m 

its undisturbed position and Tes denotes the prese: med 
mometion which depends on the mode of oscillation (j = 1; 2, 
m@ieesubScCript 7 corresponds to scattering of the regullay inei— 
Peimemwave due to the presence of the rigid body. “For the six 
iweeees Of freedom and scattering, the functions He are eeu yeie 


meoapectively, by 


a 
=> 
i 
oN 
= 
+ 
‘< 
Namen? 
= 
N 
’ 
N 
oS 
a 
iT 
N 
= 
’ 
vas 
2: 
> 
Or 
iN 
va 
= 


! [n,, Sint atiey) | + in, Cosa y || aren 


ay 
i 


pee REPRESENTATION OF THE POTENTIAL 

The solution to the boundary value problem (18) may be 
mecained by use ot a Green’s function having the physical 
macerpretation of a point wave source of unit strength. Tivese 


fources are distriputed over the suritace On tine mop cet 


- 
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according to the source strength function f sso that the 
potential at some point X, y, Z within the fluid is given by 


the surface integral 


uy cag ff fn, 0) 6 y, 2 Em, Bes (20) 
S 


where (&, n, &) represents coordinates on the surface of the 
Obyect, £f(&, n, 2) represents a distribution or weighting 
mmmetion and ds = iis (a denoece the dimensionless surface area 
element. 

The Green's function is defined, therefore, as the function 


which satisfies 


mecix, Y, 23 &, Np co) = Sx - €) Sy = TARO(z = Zz) (21) 
where 6 1s the Dirac delta function, as well as the boundary 
Memaitzions (18b, c, e). Such a function iS given by Wehausen 
and Laitone [3] as 

G(x, y, 23 & n, t) = e+ G*(x, y, z3 En, 2) (22) 


where 


Zh 


Mex, Y, 2; €, n, &) 


CO 


= fle 
? zeavs f (1+ vje¥™ cosh[u(h + n)] cosh [uth +] 5 Ginyu 
u sin u - cos u O 


(23) 
ae 2n(a’ - Vv ey cosh eit +) eos ten aii vam 
2 
ah - wv an se 


J (Cate) 


Ze 








My 7 


EEGs se ay Gy oa 8 (2 = en] 
ie onc E)? #0 Cyt ard Des, n)? +: (Zee al 
1 ae BIOS ae yz = ry2i1/2 
2- 
Vv = — = a tanh (ah) 


miger,V, indicates the Cauchy principle vale ot the intecraie 
Pumalternate series form of Green's function 1s also piven 


by Wehausen and Laitone as 


20 (v~-a*) 
mex, ¥,Z3;6,N,6) = a cosh[a(h+y)] cosh[a(h+n)]- 
avh-v ht+v 
LY , (ar) - iJ (ar) ] (Bay 
eo (uy-v") 


+4). cos [u, (h+y) ]cos[u, (htn) ]K.(u,7) 
1 uth-vh+y k k o+'k 


where Jo and 1 denote, respectively, Bessel functions of the 
first and second kind of order zero and K denotes the modi- 
fied Bessel function of the second kind of order zero. The 


quantities Wy, are the real spositive Toots oe. the equartron 


u, tan (uh) +v=0. (25) 


The solution to the boundary value problem now rests on 
the determination of the source strength function £. Taking 
the normal derivative of the potential we and applying boundary 
Somaition (18d) yields the following integral equation from 


which f is to be determined: 


OS 








1 
an ff £. (En, c)gr(x,y,258.n,5)ds = ns (5732) J. Sallie 2. ar eee 
S 


where odG/dn mav be obtained by straight forward differentiation 


@ee(22) and (24). 


fee NUMERICAL SOLUTION 

A numerical procedure can be devised by approximating, tie 
weemal Contour by a profile of a finite number OL taccus. 
Mmeeeceding in this direction, the surface of the objyectuas 
memeitioned into N area elements of size — We oles  icleKe! St.) = 
Pee takes on values j3 = 1, Z, ...N.. Since the Souree 
Semeneth function £, occurring in the integrals in (20) and 
fo) 1S a continuous, well-behaved FUNG C1ONee these mint Coma aes 


may be approximated by the following summations: 


£2... = ler (27) 
Ty ae Fal 
“hin ni? 3 (28) 
where 
mot aG : 
54 a an {f on CS Cpa ee a (29) 


(30) 


R 
il 
a om 

= 
G) 
[EN 
vad 
<< 
JH 
N 
+ 
ue 
Wi 
*< 
N 
Li, 
VV 
o. 
”n 


11 
As. 
il 


th 


Mme these expressions, wu denotes tie Demem. (aa eat eC inc an 


ni 
mraat POint on the object associated with eprener tie Neh mode 
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of o&cillatic.1, or scattering (n = 7) of the incident wave, 


The rteprations are carried out over the finite surface 
area elemen? ‘* and =n denotes the value of f at the 
cenirval pwuini-of the area element. The approximation made 


here is constptent in that as N approaches infinity, the 


approximations given in (27) and (28) approach the exact 


ROMS 

»: Mace Pegi tc es On 6, and h ere numemieal ly sevdaiudared 
(27¢*can easiely ne inverted by use of a standard computer 
subroutine to obtain the solution for f. u, is then obtained 


from (29). There are, however, certain difficulties connected 
with the evaluation of (29) and (30) which must first be con- 
_Sidered. Specifically for the special case when i = j, the 


1/Rvand o(1/h)/dn terms are singular as R approaches 0, and 


iempetondly in “the evaluation of (29) and (30) due to the singu- 


lar nature oi {i/(y tanh (uh) - a tanh (ah)] at u = a occurring 
in the infinite integral in (23) and its normal derivative. 
These difficulties are, however, circumvented by subtracting 
“poi :the singularity and carrying out its integration analytical- 
ly. These details are carried out in some, yet unpublished, 


— 


mowes fof .....-.. Garrison. 


D. HYDRODYNAM]C FORCES AND MOMENTS 

The ree and moments caused by the dynamic fluid pressure 
acting upon the immersed surface of the object may be obtained 
fem the «antesralis, 


—_ 
~ 
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Bag ies =! hee Laos a) ff Pjnyas = 4, 9) See eo (303) 


F.(t) = - (1 or a) [ J emnyas ae ee (sa 


th 


where F. denotes the i component of wave force (or moments) 


th 


and "a5 dlemIG t Co muiC mss component of forces arisine fLromneene 


th 


j component of body motion. The cocfitacrent 1-0 ac used uion 
mre case of a force (i = 1, 2, 3) whereas a is used when F 
denotes a moment (1 = 4, 5S, 6). 


Hoipunpeses Of PEECSeMEa thOlmO i NUle elUl Co mine cle slauces 


dimensionless force coefficients are defined as 


Ee | 
Cc. =e gt ee (33) 
pga”n 
and 
: eS il, 2, e 
C52 BS iG =-M5- iMG, 5 27, 2 : (oe 
pea x j by Lisg ta OF 


mmemeorresponding expressions for the moment coefficiene are 


F 


ta es eG (35) 
pga nN 
and 
J F imax 1.0) Seas ee _ os tS Cee 36) 
Ci; osa'xe 1j Me iN; | : oa, ae ( 
th 


hme complex coefficients C. eS lleycrey ice) Ke i component of 


a : ey) 
wave force (or moment) coefficient while ne denotes the it 
component of force (or moment) coefficient associated with 


the gu component of oscillation of the object. The phase 
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shift angles Fi and a5 relate the phase of the force to the 
crest of the incident wave and displacement of the body, 
respectively. The real and imaginary parts of the dimension- 
Jess force coefficients, Me and Nae are called the added 
mass and the damping coefficients, respectively. 

Using (31) and (32) an Conjunction wie (io) wanda 1 ) 
amaeche definitions (33-36), the force coefficients may be 


written in dimensionless form as 


= Jf te) 6oy,2) : eee h.(x,y,z)ds, i= 1,2,...6 (37) 
S 


and 


Ci5 - [fu;coy.2on, Goy,z)4s > yvik =Riec gee oe Oe (38) 
S 


Once pape 3 22) 1S Obtained trom (Zo), the eocttive tenecucanmp- 


G@eeaained from (37) and (38) by numerical quadrature. 


MemenaokIND'S RELATIONS AND ENERGY CHECK 

Although it may be presupposed that the numerical solution 
outlined here will converge upon increasing the number of 
Peaeeitions, modern computers are still limited by storage 
@apecity and computer time 1S expensive. In view of this, 
mes Necessary to keep the partition size as Targe, 1.2e%, 
the number of partitions 4S small, as accuracy considerations 
permit. The determination of the effect of the partition 
$size on accuracy becomes, therefore, an important considera- 


tion in order that practical limits may be established. One 
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available method is the comparison of the numerical results 
with analytical results where closed form solutions exist. 
This approach, however, is limited to a few simple shapes; 
for more general configurations no such check, of course, 
exists. : 

A method which is valid for checking the soundness of 
Mmiiesnumerical results of any Shape is the use of the so- 
called Haskind's relations as well as an energy balance, both 
of which are developed in the Appendix. The former is used 
meewerity the solution for the diffraction problem and 
mae latter the solution of the radiation problem. Using the 
asymptotic form of Green's function as given in (24) in con- 
junction with (20) the following relationship for the damping 


@eerticient 1s so obtained 


Z TT 
1 = 
me 2 Hf ¢,cenve 
oo = Vel aD ; | ee 


Beehlath+n) Je 1 2°1°°2 (8-9) G5 | 4a0 (35) 





1/2 ances =Maretane(c/e) ame biiserelaucs 


where p, = (& + z7) 
mae damping coefficient to the far field behavior of the solu- 
immo, It 15 in effect, an energy balance between the energy 
required to oscillate the body and the energy transmitted 
eeross some control surface a jJarge distance from the body. 


A relationship Somewhat similar to (39), Known @s Haskimee 


relations, may be obtained for the wave force (or moment) 
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coefficient and relates the wave produced at infinity by the 


h h 


body oscillatany fmeene ee mode to the ;° component of 


wave force as 


ae a= 
“j ~ tosh{ah) [f By Geol aS 
S 
meta (h+1) ie 1 2ereo? (Fas. (40) 


This result represents the Haskind's relation for finite depth 
as given by Seetharama Roa [4]. 

Equations (39) and (40) represent relations for the 
damping and wave force (or moment) coefficient based on the 
behavior of the far field solution. A comparison of these 
results with Ns and C, obtained from an integration of the 
Paeesure over the immersed surface, a.e€., as obtained from 
teaeand (38), provides a convenient and valuable self-check 
Om the accuracy of the numerical results. These results are 
not limited to special configurations and may be applied to 
arbitrary shapes. Equation (40) is, however, limited to 
Symmetry with respect to the x-y plane, a condition which is 


generally satisfied by practical shapes. 
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TIT. COMPUTER PROGRAM 


In order to feasibly utilize the method described in 
section D of THEORETICAL CONSIDERATIONS, a computer program 
Was caeveloped to carry out the indicated Caleulations. Although 
accuracy was a controlling factor, as it should be, the pri- 
mary requirement was to hold the computer run time to a 
reasonable value. To achieve this, every advantage was taken 
to employ the symmetry in the equations and to eliminate redun- 
dant computations by generating certain constants and matrices. 

The program consists of several subroutines, each solving 
a specific part of the problem. This allowed the development 
of the program to be simplified and the use of it was made 
easier in that only a small main program was required to 
execute the program. Subroutine GEODAT reads the input geo- 
metrical data, generates the matrix h as defined in (19), and, 
meconjunction with GCOEFF, calculates certain geometrical 
memameters. GARRIS, the primary subroutine, then generates 
the a and 8 matrices by calling either of two subroutines. 

The first (GREEN) calculates G and 3G/én based on the integral 
form given in (22) and (23) whereas the second (GREENS) makes 
evaluations on the basis of the alternate series formulations 
given in (24). For elements of the a and 8B matrices corres- 
ponding to small values of (ar) GREEN is utilized, while for 
large values of (ar) GREENS is used. With the exception of 
the diagonal elements, the majority of the elements of a and 
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B are calculated by the latter subroutine. This situation 
proved to be most fortunate because the series form converges 
rather rapidly and, consequently, requires much less computer 
time than the integral form. 

In order to check the validity of the two subroutines, 
values of G were computed. A comparison between GREEN and 
GREENS, Figs. 3, 4, 5, for constant a and decreasing r, shows 
that at small values of a, the two methods give similar 
results whereas at intermediate values of a, the GREEN solu- 
meen becomes very erratic at larger values of r. At high 
values of a, however, the two subroutines render similar 
Somutions despite the higher frequency of oscillation. This 
Stems from the fact that the amplitude of oscillation decays 
extremely fast. - 

Subroutine COMAT then determines the somrce strength 
miet1on £ by solving the complex matrix equation (27). 

Once the matrix f has been generated, the potential funceron 
mecca lculated by (28). The pressure at points on the 
immersed surface as well as the resulting forces and moments 
are then determined in subroutine FORCES. 

Eamally, the Haskind's relations and the enerey enceek are 
eimtoyed to calculate the wave force and moment cocfficients 
aude the damping coefficients in COCHCK. Here also, @ compari- 


son 1s made between these values and those found in FORCES. 
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IV S DIisCiS Sto Or, RESUIENS 


Eee VERTICALS CYR iNDER 

In order to study the accuracy of the computer program, 
mumerical results were generated for both the wave foree 
meoblem and the problem involved in the oScillationmsor othe 
object in otherwise still water. Simple geometric shapes, 
maciuding vertical circular cylinders and submerged Spheres, 
were studied and compared with closed form solutions where 
they exist. 

AS indicated in Section II, the accuracy of the iunieri-= 
Sue results is dependent on the size of the grid. Garrison 
[Ref. 5], for one, found that the fineness of the grid size, 
heemust be increased aS a increases. This resuits f£nom the 
mdiGe that the kernels of the integral equations oscillate 
rapidly as the parameter a increases. 

A Series of Computer runs was made for a vertical circu. 
mitecylindéer restting on the bottom and piercings the £ree 
Surface. Two runs consisting of cylinders with grid sizes 
he- 120 and N = 252 were made for a relative depth of h = 2:0 
mie order to determine the effect of the grid size. Another 
was made with N = 120 and h = 6.0 to determine the effect 
ee cepth. A fourth run was made with N = 240 and h = 4.0, 
Poe this time the erid was divided into © sigma ces. ) Ihe 
Dettom halt of che cylinder had N = 9O°Wath aieerenent 
aspect ratio (width/height) of 0.525 whereas the top half had 


N = 144 and an aspect ratio of 0.785.. 
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Within the range of practicality of 2) a9-50stoedas—s o50e 
the numerical results for the horizontal wave force coeffi- 
cient and the coefficient of moment about a bottom line axis 
compared favorably with the closed form solution given by 
MacCamy and Fuchs [Ref 6], Figs. 7 and 8. These results 
show that within this range of the parameter a, there is such 
a small deviation in the values obtained for the different 
grid sizes that the two curves are superimposed. As the 
frequency is increased, however, it was found that deviation 
from the closed form results increased more rapidly for 
the coarser grid than the finer one. It 1s also note- 
femeny that the magnitude of the moment coefricientmmelatay— 
to the wave force coefficient is small at low frequencies 
and increases as the frequency increases. This result is 
as expected since at small values of wavelength, wave action 
is concentrated near the free surface whereas at large 
values of wavelength, wave action is felt at deeper depths, 
thus decreasing the moment arm. 

Figures 9 and 10 represent the hydrodynamic force and 
moment coefficients associated with the cylinder oscillating 
in still water. These show that for zero frequency, the 


meee COeTricient, C approaches approximately 2m and 47 


ie? 
feecne h = 2.0 and h = 4.0 cylinder, respectively. At 

h = 6.0 a value of approximately 61 was obtained. When this 
is converted to an added mass coefficient the result is 1.0. 


Mims 1S as expected Since at zero frequiene stnemtree Siimmiec 


acts as a rigid boundary and the flow corresponds to 
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two-dimensional flow past a circular cylinder, ihe elassic a 


value for the added mass coefficient in this case is 1.0. 


Boe oPHEROID 

A series of runs for varying a was made for a half 
submerged sphere. One consisted of a grid size of 200 
points with the relative depth h = 5.0 and the others were 
made with a grid size of 264 points and h = 5.0 and 
h = 10.0. Havelock [Ref. 7] determined the added mass and 
damping coefficients for a floating, heaving sphere. His 
results were later confirmed by Kim [Ref. 8] who also 
obtained such values for floating, surging (or swaying) 
Spheroids. Alithough the results found by Havelock and Kim 
were for infinite depth, the values found in this study for 
h = 5.0 and 10.0, Figs. 11-16, compare quite favorably. 
There 1S apparently little depth effect beyond h = 5.0 
except at small values of a where it is noted the curves 
Sornesponding to h = 5.0 and h = 10.0 differ: 

Kim [Ref. 9] also determined the horizontal and vertical 
wave force coefficients as did Milgram and Halkyard [Ref. 10]. 
Comparison of these values with those calculated in this 
study, Figs. 15 and 16, show good agreement with the former, 
but some difference was noted with the latter, especially 
for the vertical force coefficient where the value of 
Milgram and Halkyard was almost three times greater at a = 2.0. 


It would appear the Milgram and Halkyard results are in error. 
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For a fully submerged spheroid there were only two 
available checks known. One is that the value of the added 
mass coefficient forja sphere in an inttprte et ude eee 
Figure 18 shows that for the numerical procedure used in 
this study a value of 0.51 was obtained for a relative depth 
6£ 5.0. The second test involves the use of Haskind's rela- 
muons Which relates @he vertical force on an oObjectutixed 
in a monochromatic wave train with the heave damping 


coefficient N55 Suen ena e 


ze «a Sia ee z : 
7D 2 2ah + sinh {2ah) &, (41) 


where N55 denotes the damping coefficient in heave. A 
comparison of the damping coefficient found in this study 
with that using (41) is shown in Table 1. Using the results 
mamene diffraction theory as the standard, the per@entage 
@eviations from it for the Haskind's relations are also 
Shown. The comparison shows that for values of a greater 
than 0.5, all of the results obtained were within 4% of the 
taken standard. For values of a less than 0.5, the differ- 
ences were considerable and this was believed to be due to 


the low value of the relative depth h. A similar compari- 


son was made for the 264 point floating sphere using h = 10.0. 


iMae results are shown in Tabiec 2. Since the present results 
Satisfy Haskind's relations, it mav be concluded the Milgrat! 
aia alkyard Ss resules are dC initelyein errer, 
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Table 1: Damping Coefficient in Heave for a Submerged Sphere 
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CC. HASKIND'S RELATIONS SA DSEVER GrsGhE Gk 

As previously indicated, important checks on the numerical 
technique and the final results can be nademby applicat tenemos 
Haskind's relations and the energy check which are developed 
m the Appendix. These checks were applied to all computer 
runs made and the percent difference, using the surface 
integration of the pressure results as the standard, is 
plotted in Figs. 19 and 20. 

In general, it was noted that dependence of the results 
Sietme grid size, N, 1S much greater for the forces obtained 
by a surface integration of the pressure when calculated by 
use of Haskind's relations. It can be seen from Figs. 19 
pa 20 that as frequency decreases, the percent deviation 
decreases. At large values of a the error becomes quite 
ieee) the effect of decreasing the grid Size or imereasine 
mae depth is also clearly shown. For the 120 point grid the 
errors are much larger than for the 240 point grid in the 
case of the h = 2.0 circular cvlinder. | Simatarly teiteme rt ome 
for the 120 point grid are much smaller for h = 6.0 than 


mor h = 2.0. 
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V SCONCES TONS 


A computer program based on linear wave theory has been 
@eveloped to calculate wave forces acline sonessupnerse 
moyects of arbitrary shape, In order to reduce thesconpura. 
meiem time, the restriction of symmetry with respect tote 
planes passing through the center of the object was 
mmicorporated. 

Oieche basis Of the results presenbed sine scetmonmins 
the following conclusions appear justified: 

im) All of €he ehecks and comparisons wenewsucecccrn™ 
mmc therefore it appears that the present method yeilds 
@ecurate results up to approximately a = 3,0 depending on 
me orid size. 

Meee line horizontal torce and moment eer. le lent Sec 
S7ses increase with the wavelength parameter a, at first, 
reach a peak and then decrease as a increases. 

See the vertical force Coefficient £o6% vay tloating 
Sphneroid starts with. a value of ti, corresponding to zero 
Size to wavelength ratio, and decreases rapidly at first 
and later slowly as a increases. 

4. Haskind's relations and the energy check give 
valuable checks on the accuracy and Show the erlect of 


miereasing the grid size or the depth. 
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APPENDIX 


HASKIND'S RELATIONS SAND ENER Gi Crime 


pee HASKIND'S RELATIONS 

The dimensionless wave force and moment components, F.(t), 
on an object of arbitrary shape may be obtained by substituting 
fom’ from (46) and (17) in (32) and simplifying. For the 


Srraction problem, this gives 
gu =A e 
F. (t) = Rela [fo i uy) xy CS € ], 2 = 1,2, 2.6 
S 


where us denotes the radiation potentials. 
The scatter potential can be eliminated from (A.1) through 
@aeceapplication of Green's reciprocal theorem, applied to G 


amd Uu. as 
J 


ff [u; (x,y, 2)0°G(x,y,258n,0)-G(X,y 256 n,0)0 U; (x,y,2)] 
R 


dxdydz = f lu, Go) GE zB) 
SetSp*S,7S 


du. 
-G(x,y,236,n, S)s(x,y,2z) Ids, 


imeethe region R, as shown in Fig. 21. Since an incompressible 
fluid was assumed, both Us and u; satisfy the Laplacian in 


in the region R, Eq. (18a), 


Bier” 
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a? OtM so ba 


a a 


(a=4) 25 


SS) 





pores -u ee = 0 (A.2) 
Sts +8, *5,, 1dn 7on 
The contributions to the surface integral in (A.2) fErom the 
surfaces Se, Sy AMG, oe WV dill S le Sai nele U4 and us Satisfy the 
iree Surface and bottom boundary conditions and the radiation 


Semadition at R=, Also, from (18d) in conjunction with (14) 


ome (15), (A.2) yields 
[| du. [| du. du, 
(uj+u) re ds = (Usa Us ae ds . (Assesa) 
S S 


Again using Green's theorem, the surface integral on the right 


hand side of (A.3) may be replaced, thus giving 


ou. ju. dU , 
[/ oop ~ “ff (US Sree on (se 
S 


S 


By using the radiation condition (18e) and since x = r €os 7 


> 


mencan be shown that 


ou. ou, 
[/ cur an Toa mn? ds = ; [f u,u, [iar(1-cos 6) 
Ss 


So 
(A.5) 


1 
= ay] ds. 


embstituting for Uy and u. from (14) and (15) and then 


Simplifying 
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etae Ju 1 
Lf (u > -u, s—) ds =.= 
o on 1 on a 
. a cosh (ah) 


jSinh (2ah) a Lar 
Gah 241 © 


aay 
"| D,(e)r-1/*[1/2-iar(1 - cos 6)] eta? 65 8 ag (A.6) 


° 


Lar «COS ic 


Inasmuch as r > ~ and e is a fluctuating quantity, 


the real part of (A.6) equals zero. Thus 


27 
- | p, (e)iar /7(1 - cos 6) Seon Ce ° dé 


© 


27 
a me, /2 mea | D; (6) (1 - cos 6) e 


° 


-ilar(l - cos 8) a6 (A.7) 


pince ar >> 1], the solution to the integral on the right 
hand side of (A.7) may be obtained by the method of stationary 


merase (Stoker [{11]). Hence 


wa ; 
nf D, (8) (1 - cos eye tar(1 - cos 8) ag 


(A. 8) 


ie 
= a Gee eae cau) D. (1) 


and after simplifying 
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du. du lye 
a O Lea 1 
[ CE arma u, ae ds = Lia) See 
2 cosh’ (ah) 


j sinh (Zah) 
2a 


+ h] e!T D. (1) oS) 


From (A.1),. (A.3)5 (A.4) and (AS), simpliticaticimew e- 


lef 
.(2Ta) sinh(2ah) 


Sill 
Meee) = Rel-i at Gre ti (A.10) 
‘i cosh” (ah) ae 


and taking the modulus on both sides gives 


1/2 
By (27a) r Sinh (2ah) 


i | Cet boa lie 255.58 (A.11) 
: cosh“(ah) oe : 


Im actuality, (A.11) represents six difference canatian. 


which are known as "Haskind's relations."' They relate the 
jth component wave force (or moment) coefficient to the waves 
th 


pmeauced at infinity by the object oscillating in the i moa ce’. 


bem ENERGY CHECK 


en component of the dynamic 


Eiilaciom (5.) represents Ene m1 
Herce {Or moment) arising due to the pen mode of oscillation . 
Oemeene Object. In the case of the radiation problem, they 
may each by expressed as the sum of two components, one in 


phase with the acceleration and the other in phase with the 


velocity of the body in the form 


Fy, (t) = - My,ax,(t) - Nj ;aX;(t) , i,j = 1, 2,...6 (A.12) 


where M represents the added mass (or added moment of 


inertia) and is defined, when compared to (34), (36) and (37), 


OZ 





as 


—_ 
— 
if 


- 3 1 
-oa rel ff uj(x.y,2)h Gy,2)ds] Mey eee (A.13a) 
S 


mem a force and for a moment as 


4 1. “= 4555 
-pa ret f [uj Ccy.2)hy @oy,2)4s}, Me DG. (A, 1359 
S 


ij 
and Ni denotes the damping coefficients (in relation to 
energy transport) and is defined, when compared to (34), (36) 


ame (37), as 


: a 1, =o See 
ij = -00a nt f f aj Ccy,2)hy Gy, 2)d5], fe 1. 2 6 (A.14a) 
S 
home a force and 
7 _A 1 = eo 
“a -o0a int ff a; cx.y,z)h; (x,y, 2)4s], HS dn Pooch (A.14b) 
S 
hone a moment. 
th 


ine dimensional £orce (or moment) component, in the 1 
aire ction, on the oscillating body due to oscillations in the 
a mode may then be obtained from (A.12) by equating j to 1. 
pemmee the forces and moments exerted by the body on the fluid 
are equal and opposite to these, the energy transmitted by 


the object to the surrounding fluid during one period is given 


by 
T 


E. = | Gh or 2) Be, Ce) X. (t) dt (ALS) 


mitre a 15 used tor a force (1 = 1,92 S)eandeene factor me ton 


aemoment (i = 4, 5S, 6). After carrying out the integration, 


os 





(A.15) may be rewritten as 


Eee = Tce Ne 
1 1 


' i, <2 eo (A.16) 


i? 
where NG denotes the dimensionless damping coefficient which 


moeaefined for a force as 


N.. 1 = 3 
g 0) b > 
7 ’ a se 
1) 00a J 3 es a6 


< ij i= 4, S, 6 
ij ~4 
J 00a 3 





Jen he eee Or 


\leg n; is the amplitude of the outgoing progressive waves 


aun mode, then 


aioe aue to oscillation of the body in the i 
fees amplitude can be related to the asymptotic velocity 
mocential, Us; by using the dynamic tree surtace bouldiaa, 


Gendition and 


nn (GegGh) = aX? tanh(kh) a Ju, (r,6,0) | ; (Aca) 


The total energy transmitted over one period can now be 





written 
: - -4 4 
E. bt Ag [1 + Z2kh = aaG ee r 
exe sinh(2kh) §& 
Zi 
| Ju; (r,8,0)|%ae t= 1 eee oe (A.18) 


° 


Now, by comparing (A.16) and (A.18), the damping coeffi- 


cients may be written in terms of the asymptotic velocity 
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Mote mt wai le U;, as 


~ bri 2ah 
roe 2 (onan) anne 


ri 


f ju, GOs 0)itede |) i= tau emer CoN) 
° 


C. NUMERICAL EVALUATION OF HASKIND'S RELATIONS AND ENERGY CHECK 
The right hand sides of (A.11) and (A.19), involving the 
asymptotic velocity potentials, are solved by the use of the 
@esceribution functions f., obtained from the numerical solu- 
tion of the radiation problem, in conjunction with the series 
morm Of Green's function (24). 
For eed and y = 0, considerable simplification is possible 


in the asymptotic form of G and may be written as 


Qn (a--v*) 
ee s2565",6) = ~ = —>— cosh(ah) coshfa(n-h)] 
a h-v htv 
eee 
poe Fe (A.20) 


Tar 
1 


Thus the asympecotieG potentials Us; 1S given by 


Py tee Z 1/2 (v*-a*)cosh (ah) 
u,(r, 0, ) =A os 5 aay ? 
a h-v htv 
Nar) 
If £.(é, mW. Cc) GOsn [a(n+h)]r, e 1 4°ds 2 Gea) 
S 
bimec Tt ce large, r, may be replaced by r and moved outside 

Mme integral, Further, as shown in Fig. 22, £o0r large F 


Be Ey cos(B-86) 
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(2 ae 


x) 


LOATdO@ Jno aie Ny a 


°7Z oin3sty 
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IZ 


where Pp, = ae + 7°) eee! u, may be rewritten as 
Sly Le oe ly A 7s eee 7 
u(r, oT yee? as (ve-dn cosh (al) o oi(ar 70) 
a -h-v htv 
[ f.(E, n, )cosh[a(nth)] pues) ye (Anz) 
S 


substituting (A.22) into (A.19) and simplifying gives the 


result 


1 a=" 


Tl 
ia wy a ee ; SJ ae To esvieos mn fant be 


(e714 cos (B-8) - e 1ap,cos (B+8)) a. 2 de. (A.23) 


Once u,(T, 6, 0) is known, D. (7) may be obtained by using 
the radiation condition (18e). Thus 


1/2 1ar 
D. (7) = u;(r,7,0)r } (A.24) 


ameearter simplifying by means of trigonometric identifies, 


(A.11) becomes 


a orcs | £;(E,v,t)cosh[a(nth) Je *#71°°S "Mas, (A.25) 


ieZs 
wiltere 1 = 1, 2,°°6 since for an axisymmetric body, the only 
muses Of interest are those pertaining to i= 1, 2, and 6. 


A comparison of the results obtained from (A.23) and (A.25) 
with those calculated by integration of the pressure over the 
submerged surface, Eqs. (37) and (38), are a valuable seltr- 


check on the accuracy of the numerical results. 
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